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Abstract 

We analyze the effective triplet interactions between the centers of star poly- 
mers in a good solvent. Using an analytical short distance expansion inspired 
by scaling theory, we deduce that the triplet part of the three-star force is 
attractive but only 11% of the pairwise part even for a close approach of three 
star polymers. We have also performed extensive computer simulations for 
different arm numbers / to extract the effective triplet force. The simulation 
data show good correspondence with the theoretical predictions. Our results 
justify the effective pair potential picture even beyond the star polymer over- 
lap concentration. 
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I. INTRODUCTION 



Star polymers i.e., structures of / linear polymer chains that are chemically linked 
with one end to a common core, have found recent interest as very soft colloidal particles 
0-^. As the number / of chains increases, they interpolate between linear polymers and 
polymeric micelles For large /, the effective repulsion between the cores of different 

polymer stars becomes strong enough to allow for crystalline ordering in a concentrated star 
polymer solution. While such a behavior was already predicted by early scaling arguments 
only recently corresponding experiments have become feasible with sufficiently dense 
star solutions. The crystallization transition occurs roughly at the overlap concentration c* 
which is the number density of stars where their coronae start to touch experiencing the 
mutual repulsion. It is defined as c* = l/{2RgY where Rg, the radius of gyration, is the root 
mean square distance of the monomers from the center of mass of a single star. In addition, 
theory and computer simulation have refined the original estimate for the number of chains 
/ necessary for a freezing transition from / ~ 100 to / ~ 34 ||^ and predicted a rich 
phase diagram including stable anisotropic and diamond solid structures at high densities 
and high arm numbers. These results were derived using an effective pair potential between 
stars with a logarithmic short distance behavior derived from scaling theory. 

In general, while the pair interactions are the central focus and the typical input of any 
many-body theory, much less is known about triplet and higher-order many body interac- 



tions. For rare gases, the Axilrod- Teller triplet interaction [|IU| has been found to become 



relevant in order to describe high-precision measurements of the structure factor For 



charged colloids, the effective triplet forces are generated by nonlinear counterion screening. 



This was investigated recently by theory and simulations |]I2|. For star polymer solutions 
in a good solvent such studies are missing. In all three cases, the effective triplet forces 
originate from formally integrating out microscopic degrees of freedom. For rare gases, 
these are the fluctuations of the outer-shell electrons while for charged colloids the classical 
counterions play the role of additional microscopic degrees of freedom. For star polymers, 
on the other hand, one is interested in an effective interaction between the star centers by 
integrating out the monomer degrees of freedom |T^. Usually one starts from an effective 
pair potential which is valid for large particle separation. The range of this effective pair 
potential involves a certain length scale £ which is the decay length of the van-der-Waals at- 
traction, the Debye-Hiickel screening length or the diameter of gyration 2i?g, for rare gases, 
charged colloids, and star polymers, respectively. Triplet forces, i.e. three star forces, not 
forces between monomers, become relevant with respect to the pairwise forces if the typical 
separations between the particles are smaller than this typical length scale £. This implies a 
triple overlap of particle coronae drawn as spheres of diameter £ around the particle centers. 
The triple overlap volume is an estimate for the magnitude of the triplet forces. Hence 
a three-particle configuration on an equilateral triangle is the configuration where triplet 
effects should be most pronounced. 

The aim of the present paper is to quantify the influence of triplet interactions for star 
polymer solutions in a good solvent using both analytical theory and computer simulation. 
In doing so, we consider a set-up of three star polymers whose centers are on an equilateral 
triangle. We found that the triplet part is attractive but its relative contribution is small 
(11%) with respect to the repulsive pairwise part. This relative correction is universal, i.e., 
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it is independent of the particle separation and of the arm number. It even persists for a 
colhnear configuration of three star polymers where the absolute correction is smaller than 
in the triangular situation for the same star-star distance. Consequently, the validity of the 
effective pair potential model is justified even at densities above the overlap concentration. 
In particular, our result gives evidence that the anisotropic and diamond solids predicted 
by the pair theory are indeed realizable in actual samples of concentrated star polymer 
solutions. 

Our paper is organized as follows: in section II we apply scaling theory to extract the 
triplet forces both for small and for large arm numbers. In section III we briefly describe 
our Molecular Dynamics (MD) simulation scheme and present results in section IV. Com- 
paring these to the theoretical predictions, we find good agreement. Section V is devoted to 
concluding remarks and to an outlook. 



II. SCALING THEORY OF TRIPLET FORCES BETWEEN STAR POLYMERS 



A. Scaling of single stars 

The scaling theory of polymers was significantly advanced by de Gennes' observation that 
the n-component spin model of magnetic systems is applicable to polymers in the formal 
n = limit [|1^]. This opened the way to apply renormalization group (RG) theory to explain 



the scaling properties of polymer solutions that have been the subject of experimental and 
theoretical investigations since the pioneering works in this field [|I5|. Many details of the 
behavior of polymer solutions may be derived using the RG analysis [|16|. Here, we use only 
the more basic results of power law scaling: the radius of gyration Rg{N) of a polymer chain 
and the partition function Z{N) are found to obey the power laws: 

Rg{N) N" and Z{N)^z^N^-\ (1) 

The fugacity z measures the mean number of possibilities to add one monomer to the chain. 
It is microscopic in nature and will depend on the details of the model or experimental 
system. The two exponents u and 7 on the contrary are the n = limits of the correlation 
length exponent z/(n) and the susceptibility exponent 7(n) of the n component model and 
are universal to all polymer systems in a good solvent, i.e., excluding high concentration of 
polymers or systems in which the polymers are collapsed or are near the collapse transition. 
For any such system the exponents of any other power law for linear polymers may be 
expressed by these two exponents in terms of scaling relations. 

It has been shown that the n component spin model may be extended by insertions of so 
called composite spin operators that allow to describe polymer networks and in particular 
star polymers in the n = limit |1T7|- [19|] . A family of additional exponents 7/ governs the 
scaling of the partition function Zf{N) of a polymer star of / chains each with N monomers: 

Zf{N) ~ z^N^f-\ (2) 

Again the exponents of any other power law for more general polymer networks are given by 
scaling relations in terms of 7/ and u. Here, we substitute another family of exponents ?7/ to 
replace 7^ — 1 = z/(?7j — /?72)- The first two members rji = and 172 = (1 ~ 7)/^^ are defined 
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by the requirement that the / = 1 star and the / = 2 star are just hnear chains with scahng 
exponents 71 = 72 = 7- The values of these exponents are known from renormahzation 
group analysis (RG) and Monte Carlo (MC) simulations . Several equivalent approaches 
have been elaborated to evaluate the renormalized perturbation theory. Early first order 
perturbative RG results were given in ref 



expansion in the parameter e 
the rjf reads 



2T[] . Here, we explicitly present the result of an 
A — d where d is the space dimension. The e-expansion for 



-!/(/ 



~U^f- 25) + ^ [(28/ - 89)C(3) + 8/^ - 49/ + ^] } + 0{e') (3) 



with the Riemann ^-function. Note that this series is asymptotic in nature and to evaluate 
it for e = 1 it is necessary to apply resummation. An alternative expansion for the star 
exponents makes use of an RG approach at fixed dimension d = ?> proposed by Parisi p2 . 
This expansion has been worked out in refs. |]23|- p5|] . The corresponding expressions are 
lengthy and not presented here. In Table 1, in the first two lines we have calculated the 
resummation for the series in Eq. (|^) as well as for the expansion at fixed dimension. 
The resummation procedure that we apply combines a Borel transform with a conformal 
mapping using all information on the asymptotic behavior of the perturbation expansion 

. Results for / < 9 have been given before in refs. 
whereas we have added here the calculation of values for / = 10, 12, 15. The 
deviation between the two approaches measures the error of the method. For large / the 
leading coefficient of the fcth order term in Eq. is multiplied by f^^^. This is due to 
combinatorial reasons and occurs also for the alternative approach. It limits the use of the 
series to low values of /. 

Another possibility to estimate the values of the star polymer scaling exponents 7/ is to 
consider the limiting case of many arm star polymers. For large / each chain of the star is 
restricted approximately to a cone of solid angle Qf = Air/f. In this cone approximation 
one finds for large / 



of the corresponding spin model p6 



7/ 



(4) 



B. Two star polymers 

Let us now turn to the effective interaction between the cores of two star polymers at 
small distances r that are small on the scale of the size Rg of the stars. Let us for the moment 
consider a more general case of two star polymers with /i and /2 arms respectively. The 
cores of the two stars are at a distance r from each other. We assume all chains involved to be 
of the same length. The power law for the partition sum ^j-Jjj ^^^^ polymers may 

then be derived from a short distance expansion. This expansion is originally established in 
the field theoretic formulation of the n component spin model. While we do not intend to 
give any details of these considerations here, applications to polymer theory may be found 
in refs. [^,^. The relevant result on the other hand is simple enough: the partition sum 
of the two stars Z^^j^{N, r) at small distance r factorizes into a function Cf^f^{r) of r alone 
and the partition function Zf^^f,^{H) of the star with /i + /2 arms that is formed when the 
cores of the two stars coincide. 
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Zu.{N,r)^Cf,f,{r)Zs,^f,{N) (5) 
For the function Cf^f^{r) one may show that power law scahng for small r holds in the form 

CMAr)-r''M2. (6) 

with the contact exponent 0}/j2 ■ ^'^'^ scaling relation for this power law we we change 
the length scale in (|^) in an invariant way by r — > Ar and — > X^^'^N. The scaling of the 
partition function -Zj^jj may be shown to factorize into the contributions for the two stars. 
This transforms (|]) to 

^-1/^7/1 -1) ^-1/^7/2 -i)2(2^^(^i/^']V,Ar) ~ X'^M2Cf,f,{Xr)X-^/''^^fi+f2-^'>Zf^+f^{X^/''N) . (7) 
Collecting powers of A provides the scaling relation 

^0/?/. = (7/1 - 1) + (7/. - 1) - (7/1+/. - 1) , 

®/?/2 = + Vf2 - Vh+f2 ■ (8) 



of arms /i = /2 = /• The mean force Fjj{r) between the two star polymers at short distance 

r is then easily derived from the effective potential V'^^{r) = —kBT\og[Z^^j(r)/(Zf)'^] with 
k^T denoting the thermal energy. For the force this results in 

-1-F;j)(r) = ^. (9) 



The cone approximation for the contact exponents [^] G^^^ may be matched to the known 
values for / = 1, 2 (see table 1), fixing the otherwise unknown prefactor. Assuming that the 
behavior of the 0^^ may be described by the cone approximation for all / one finds: 

^f/^(0-^V-- 

This matching in turn suggests an approximate value for the rjf exponents, 

^/--^(2^/^-2rr/^ (11) 

Note on the other hand that this approximation is inconsistent with the exact result rji = 0. 
However, the approximation works well for Qjj in the range / = 1, . . . , 6 were we have 
calculated the corresponding values from the perturbation theory results as well as according 
to the cone approximation. Our results, displayed in the second part of table 1, show good 
correspondence of the cone approximation with the resummation values. 



5 



C. Three stars 



We now use the idea of the short distance expansion once more to derive the triplet 
interaction of three star polymers at close distance. We consider a symmetric situation 
in which the three cores of the polymer stars are located on the corners of an equilateral 
triangle (see Fig. 1). The distance between the cores is r while their distance to the center 
of the triangle is R. We assume that the radius of gyration Rg of the star polymers is much 
larger than their mutual distance Rg^ r. 

To make the argument more transparent we first consider the slightly more general case 
of three stars with /i, /2 and /s arms respectively. Shrinking the outer radius R of the 
triangle on which the cores are located, the partition function of this configuration of three 
stars will scale with R according to 



Zf,j,f,{R)-R^hf^f^ (12) 
'hhh 



Q/1/2/3 = ^/i + + % - ^/i+/2+/3 • (13) 



Now, the scaling exponent ?7/i+/2+/3 of the star that results by collapsing the cores of the 
three stars at one point has to be taken into account as follows from an argument analogous 
to the above consideration for two stars. 

Let us specify the result for the symmetric situation of three equivalent stars /i = 
/2 = /s = /• Furthermore we assume that the large / approximation (|TT]) is valid for the 
exponents rjf. Then the three star contact exponent may be written as 

An effective potential of the system of the three stars at small distance R from the center 
may then be defined by 

V^ffiR) = -kBTQfjf HR/Rg) . (15) 

We now derive the corresponding three body force underlying this effective potential. Note 
that the absolute value of the force is the same for all three stars. The relation of the 
potential to the force on the core of one star is then 



vjf;''{R + dR) - vfffiR) = Y^F,-dR, = 3Fj%{R)dR . (16) 

1=1 



The final result for the total force on each of the stars that includes any three body forces 
is therefore 

Fj%iR) = -k^Tefjj/iSR) . (17) 

If one starts instead from a sum of two body forces, then one star experiences the sum 
of the two forces calculated for the star-star interaction. With the given geometry of the 
equilateral triangle this is easily calculated to give 

Fi%{r) = IruOff/ru + hsOff/rnl = -k^TQ^j/R. (18) 
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Here, r = ri2 = = R\/3 denote the distance between two of the stars, while the fij are 
the unit vectors along the edges of the triangle (see Fig.l). The relative deviation from the 
pair potential picture is then given by 



(3) (r) - F^^) (r) B^^) - 38^^^ 



^fff ^fffy) '^^ff 



Using the cone approximation for the contact exponent we finally obtain for the relative 
deviation caused by triplet forces alone 

AF 33/2 - 3 

^-^--O-ll- (20) 

This result is independent of the number of arms and valid in the full region that is described 
by the logarithmic potential. In table 1 we have calculated the exponents as derived from 
the perturbation expansion of polymer field theory pU| , P^p5| checking the relation eq. i^^- 



Taking into account the error that may be estimated from the difference of the results 
obtained by the two complementary approaches, the results are in good agreement with 
the cone approximation even for low / values. The fair coincidence is rather surprising as 
additional numerical errors might be introduced by the calculation of the contact exponents 
from the original star exponents. It confirms our estimate of the relative deviation caused 
by triplet forces to be of the order of not more than 11% for all analytic approaches we have 
followed here. Let us note that the analogous calculation for a symmetric linear configuration 
of three stars yields the same relative deviation eq. (|20|) . The absolute triplet forces for the 
linear configuration are smaller by a factor -\/3/2 than for the triangular configuration with 
the same star-star distance. 



III. COMPUTER SIMULATION METHOD 

Molecular dynamics (MD) simulations were performed using exactly the model that 
three of the present authors devised to test the effective pair potential ^ and had been 
originally proposed to study single star polymers p2| , |33[] . In this model the configuration of 



star polymer i = 1, 2, 3 is given by the coordinates r^'-'-* of the N monomers m = 1, . . . , of 
the / chains j = 1, . . . , / and the position of its core The main features of this model are 
the following: (1) A purely repulsive truncated Lennard- Jones like potential acts between 
all monomers m = 0, . . . , on all chains. (2) An attractive FENE-potential |]32| , |33| that 



preserves the chain connectivity and acts only between consecutive monomers m, m + 1 along 
each chain. (3) These potentials have to be slightly modified for the interaction between the 
first monomer m = 1 and the core m = of the star to allow the core to have a radius 
that is sufficiently large to place / monomers in its vicinity. 

The three cores of the stars were placed at the corners of an equilateral triangle, see 
again Figure 1 where also the core radius Rd is shown. A typical snapshot of the three star 
simulation is displayed in Figure 2 for a functionality of / = 5 and = 100 monomers per 
chain. The force on the star core was averaged during the MD simulation for a number of 
edge lengths r of the triangle varying in the range between the diameter of the two cores 
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2Rd and the diameter of gyration 2Rg of a single star polymer. We have produced data for 
/ = 3, 5, 10, 18, 30. For the smaller functionalities (/ = 3, 5, 10) the number of monomers 
per chain was = 100 while for / = 10, 18, 30 a number = 50 was chosen. Note that the 
total system comprises between 900 — 4500 mutually interacting particles. As equilibration 
is slow and the statistical average converges slowly, the simulation becomes increasingly 
time-consuming beyond such system sizes. As for reference data, we have also produced 
data for a two stars situation according to the calculations in Ref. [Q. 

IV. RESULTS 

Results of the computer simulation are compared to the theory in Figures 3a and 3b. 
The reduced averaged force on a single star is shown versus the reduced triangle length for 
different arm numbers. As a reference case, also the corresponding results in a pair potential 
picture are shown, both within theory and simulation. For technical reasons we kept a small 
core radius Rd in the simulation, which is roughly 10% of the radius of gyration of the whole 
star. In the theory, on the other hand, the core size was zero. Hence, to compare properly 
0, a shift r — 2Rd has to be performed. 

As expected, in both theory and simulation, the triplet forces become relevant only within 
the coronae. A comparison with pure pairwise forces leads to the first important observation 
that the triplet force is smaller, i.e. the pure triplet contribution is attractive. (Note that one 
has to multiply the pure two-star force by a factor of y/S for simple geometrical reasons.) 
The relative magnitude of the triplet term, however, is small. A quantitative comparison 
with theory and simulation leads to good overall agreement. The triplet contribution itself, 
however, is subjected to larger statistical errors of the simulation. Hence we resorted to 
a different strategy to check the theory by plotting the inverse force versus distance. If 
the theory is correct the simulation data should fall on a straight line both for the pure 
pairwise and the full triplet case. The slope should then give the theoretical prefactor of the 
logarithmic potential. The advantage of this consideration is that the slope bears a smaller 
statistical error as more data points are included. Such a comparison is shown in Figure 4 
for / = 10. The first consequence is that the simulation data indeed fall on a straight line 
confirming the theory. In fact this is true for all other parameter combinations considered 
in the simulations. The slope is higher for the triplet and lower for the pair case, both in 
theory and simulation. The actual values in Figure 4 are in the same order of magnitude 
but a bit different. 

In order to check this in more detail, we have extracted the slope for all simulation data. 
The result is summarized in Figure 5 where the relative differences of the slopes between the 
pair and triplet cases are plotted versus the arm number /. The theory predicts a constant 
value of 0.11, see Eq. (^). The simulation data scatter a lot in the range between 0.05 and 
0.15 due to the large statistical error but the theoretical value falls reasonably within the 
data. Consequently, the triplet contributions are found to be attractive and small even for 
nearly touching cores where the triplet overlap of the coronae is substantial. 
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V. CONCLUSIONS 



In conclusion, we have calculated, by theory and computer simulations, the triplet in- 
teraction between star polymer centers in a good solvent positioned on the corners of an 
equilateral triangle. The triplet part was found to be attractive but only about 11% of the 
pairwise repulsion. Our calculations justify earlier investigations where the pair potential 
framework was used even slightly above the star overlap concentration. 

We finish with a couple of remarks: First, the scaling theory can also be performed for 
any triplet configurations beyond the equilateral triangle studied in this paper. Second, 
arbitrary higher-order many body forces can be investigated assuming a cluster of M stars. 
Such a calculation is given in Appendix A. As a result, the deviations from the pair potential 
picture increase with the number M and even diverge for M —>■ oo. This implies that the 
pair potential picture breaks down for very high concentrations. This is expected as for 
high concentration a star polymer solution is mainly a semi-dilute solution of linear chains 
where it is irrelevant at which center they are attached to . As far as further simulational 
work is concerned, there are many open problems left. Apart from the investigation for 
arbitrary triplet configurations and their extensions to an arbitrary number of stars, the 
most challenging problem is a full "ab initio" simulation of many stars including many-body 
forces from the very beginning. This is in analogy to Car-Parrinello simulations which 
were also applied to colloidal suspensions [^. A first attempt has been done [0, but 



certainly more work is needed here. Another (a bit less demanding) task is to study stars 
on a periodic solid lattice with periodic boundary conditions and extract the many body 
interactions from there. 

It would be interesting to study the relevance of triplet forces for star polymers in a poor 
solvent near the B-point [^. It can, however, be expected that the triplet forces here are 
even less important than for a good solvent as the effective interaction becomes stiffer in a 
poor solvent. Furthermore, the effect of polydispersity in the arm number which has been 
briefly touched in our scaling theory treatment should be extended since this is important 
to describe real experimental samples. 
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APPENDIX A: HIGHER ORDER FORCES BETWEEN STAR POLYMERS 

Here, we derive for the general case of M simultaneously interacting star polymers with 
/ arms the effective Mth order force. Generalizing the equilateral triangle geometry, we 
study the situation where the M cores of the stars are evenly distributed on a sphere with 
radius R. In particular, the cores of the stars may be located at the corners of a regular 
polyhedron. Then the non-radial forces on each star polymer cancel. The latter condition 
may be fulfilled approximately also for large numbers M for which a regular polyhedron 
does not exist. 
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We first calculate the force on one star by the sum of M — 1 pairwise forces effected by 
the other stars. For the pairwise force @ that acts according to a 1/ r-law it is easy to verify 
that the radial component of the force between any two points on the sphere is 0jj/(2i?) 
independent of their relative position. With this simplification the total (radial) force on 
one star is 



kBT ^'^ 2 R 

(2) 

Here, -^m / denotes the sum of pairwise forces on one of the M stars each with / arms. In 



the case M = 3 this is the result of eq. (p^S]). 

The total Mth order force F^f^j between M star polymers with / arms brought close 
together may again be derived from a short distance expansion resulting in the scaling 
relation 

eiyj = M.rjf- rjM.f . (A.2) 
The force on one star is then found in the same way as for three stars as 



M-R 

The leading contributions for large numbers of stars M in the two cases differ even in the 
power of M. While the first is linear in M the latter grows only with the square root of M. 
In the large-/ and large-M approximations this reads : 

' ^fS!, « (A.4) 



keT ~ 18 2R 

Note that for large M the factors M and M^/^ in these two approaches are not a result of 
the large / approximation but are of combinatorial and geometrical origin. This shows that 
for large M the sum of pairwise forces largely overestimates the force on one star. 
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TABLES 



f 


1 


2 


3 


4 


5 


6 


8 


9 


10 


12 


15 


arif 
b 






-0.28 
-0.28 


-0.75 
-0.76 


-1.36 
-1.38 


-2.07 
-2.14 


-2.88 
-3.01 


-4.71 
-5.06 


-5.72 
-6.22 


-6.80 
-7.48 


-9.12 
-10.23 


-12.98 
-14.93 


a efj 

h 
c 


0.28 
0.28 
0.28 


0.80 
0.82 
0.79 


1.38 
1.49 
1.44 


1.99 
2.30 
2.22 


2.66 
3.20 
3.11 


3.36 
4.21 
4.08 












(3) 
^ ^fff 

b 
c 


0.75 

0.76 
0.74 


2.04 

2.17 
2.08 


3.47 

3.94 
3.83 


5.04 

6.09 
5.89 


6.77 
8.51 
10.82 














a AF/F 

b 

c 


-0.11 
-0.09 
-0.11 


-0.15 
-0.12 
-0.11 


-0.16 
-0.12 
-0.11 


-0.16 
-0.12 
-0.11 


-0.15 
-0.11 
-0.11 















TABLE I. Calculation of the exponents that govern the pair and triplet interactions. The 
labels (a) and (b) stand for the two complementary renormalization group approaches (expansion 
in e = 4 — d and massive renormalization at d = 3) used to calculate the exponents rjf. The 
difference of the two results may be taken as an estimation of the error of the method. Label (c) 
stands for the cone approximation result with matching to / = 1, 2 as explained in the text. 
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APPENDIX: FIGURE CAPTIONS 



FIG. 1. Three star polymers at mutual distance r. The cores of the stars (with radius 
Rd) are located at the corners of an equilateral triangle. The distance from the center is R. 
The mean radius of gyration of a single star is Rg. 

FIG. 2. Snapshot of the simulation of three stars with / = 5 arms each with = 100 
monomers. The cores are located at the corners of the equilateral triangle that is depicted 
in the center. The monomers that belong to the same star are represented by balls of the 
same color: either black, dark gray, or light gray. 

FIG. 3a. Comparison of the force F measured in the three star MD with that calculated 
from a corresponding two star MD simulation for / = 3 and / = 10 with N = 100. Also 
the results predicted by the theory are plotted as a continuous line (only pair forces) and a 
broken line (including triplet forces). 

FIG. 3b. Same as Fig. 3a but for / = 18 and / = 30 with N = 50. 

FIG. 4. Comparison of the inverse force 1/F measured in the three star MD with that 
calculated from a corresponding two star MD simulation for / = 10 with N = 50. The linear 
fits for the pair forces (small dashed hne) and the full three body force (dash-dotted hne) 
are shown together with the respective results predicted by the theory which are depicted 
by a continuous line (only pair forces) and a broken line (including triplet forces) . 

FIG. 5. The slopes of the linear fits to the data as shown in Fig. 4 were extracted from 
the simulation data for / = 3, 5, 10, 18, 30 and N — 50, 100 to calculate the relative deviation 
AF/Pi'ff induced by the triplet forces. The hne at 0.11 corresponds to the analytic result. 
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Fig.l, von Ferber et al. 
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Fig.3a, von Ferber et al. 




Fig.3b, von Ferber et al. 




Fig. 4, von Ferber et al. 
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